1 



(N 
> 

o 



> 
in 

o 



(N 



X 



Charmed Meson Spectroscopy and Matrix Elements with an 
0(a)-Improved Clover Fermion Action 

UKQCD Collaboration 
presented by David G. Richards 

Department of Physics, University of Edinburgh 
Mayfield Road, Edinburgh EH9 3JZ, Scotland 

We present preliminary results for the spectrum and decay matrix elements for heavy-light and heavy-heavy 
mesons, obtained on the 64-node Meiko Computing Surface at the University of Edinburgh. 

Quark propagators are computed with an 0(a)-improved fermion action on 24^ x 48 lattices at /3 = 6.2, using 
three values of the quark mass up to around the strange quark mass, and four values of the quark mass in 
the region of the charm quark mass. We compare results for the hyperfine splitting in charmonium with those 
obtained using the conventional Wilson fermion action and find that the splitting is 1.83(15) times larger with 
the improved action. Our measurements of fs indicate non-scaling corrections of the order of 20% to the Heavy 
Quark Effective Theory expectation. A comparison is made with results obtained on 16'^ x 48 lattices at /3 = 6.0. 



1. INTRODUCTION 

This contribution represents the talk "Latest 
Results from the UKQCD Collaboration" pre- 
sented at Lattice 92. 

As the masses of hadrons in lattice gauge simu- 
lations approach the inverse lattice spacing, a~^, 
quantities are expected to become increasingly 
sensitive to discretisation errors. Thus it is in this 
regime, corresponding in present simulations to 
quarks with masses of the charm and above, that 
we expect differences between an 0(a)-improved 
fermion action and the standard Wilson fermion 
action to be most apparent. 

In this talk I shall describe the recent studies 
by the UKQCD Collaboration of mesons contain- 
ing propagating charm quarks. The talk begins 
with a brief description of the computational de- 
tails. Then a comparison is made of the hyper- 
fine splitting in charmonium between the stan- 
dard Wilson action, and the 0(a)-improved clover 
action For the clover action, this investiga- 
tion is extended to the D and B systems. Finally, 
measurements are made of the pseudoscalar decay 
constant, fp, for heavy-light systems. By deter- 
mining the behaviour of jp^Mp as a function of 
the pseudoscalar mass, Mp, an estimate is made 
of the magnitude of the non-scaling terms in the 



Heavy Quark Effective Theory (HQET) in the re- 
gion of the i?-meson mass. 

2. COMPUTATIONAL DETAILS 

Clover quark propagators were generated on 
quenched 24'^ x 48 lattices at /3 = 6.2, using 
three values of the quark mass up to the re- 
gion of the strange quark mass, with hopping pa- 
rameters Ki = 0.14144,0.14226 and 0.14262 [||. 
In addition, propagators were generated at four 
values of the quark mass in the region of the 
charm quark mass, with hopping parameters 
Kh = 0.133,0.129,0.125 and 0.121; the charm 
quark mass corresponds to K/i ~ 0.129. The 
preliminary results presented in this write-up are 
based on an analysis of 33 configurations, the re- 
sults in the talk being based on an analysis of 21 
configurations. 

The light quark propagators were generated 
using local sources and sinks (LL). The charm 
quark propagators were computed from smeared 
sources, to both local sinks (LS) and smeared 
sinks (SS), using a Jacobi-smearing algorithm ||] 
with a smearing radius r ~ 5. 

The comparison of the hyperfine splitting 
in charmonium between the Wilson and clover 
fermion actions was obtained on a common sam- 
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pie of 18 configurations, using local sources and 
sinks [Q . All fits to the correlators are performed 
taking into account the correlations between dif- 
ferent time slices, with the error bars obtained 
using a bootstrap procedure. 

3. HYPERFINE SPLITTING 

Quenched calculations with the standard Wil- 
son fcrmion action typically find values of vector- 
pseudoscalar mass splittings, to/j. — m^i and 
rrijj^ — niji^, which are much too small, indeed 
typically by a factor of 2 and 4 respectively at 
/3 = 6.2 1^,^. The clover fermion action is related 
to the standard Wilson fermion action through 
the addition of a term. 
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The hyperfine splitting is expected to be sensitive 
to the magnetic moment form of this additional 
term. 

To compare the values for the hyperfine split- 
ting in charmonium, we chose k = 0.135 for the 
Wilson action, and k — 0.129 for the clover ac- 
tion, yielding approximately the same value of 
the pseudoscalar mass, Wr/ci fo^' the two actions. 
The measured values of the vector and pseu- 
doscalar masses are shown in Table 0. Setting the 
scale from the string tension gives = 2.73(5) 
GeV and we find the values of the masses 
quoted in the table to be within a few percent 
of their physical values. However, the mass split- 
ting quoted in the third column is very different 
for the two actions, and we find that 



(m 



J/i: 



clover 



= 1.83 



13 
15 



(2) 



Using the value of a^^ quoted above, we obtain 
mj/^-m^^ = 281: i MeV (Wilson) 
mj/^ — m^^ — 521 ^ MeV (clover). (3) 

to be compared to the experimental value of 
117(2) MeV. The corresponding values found by 
El-Khadra et al. 0, using a mean field calcula- 
tion of the coefficient of the clover term in eq. |l|, 
are: 51(3) MeV at /3 = 5.7, 62(4) MeV at /3 = 5.9 



and 68(5) MeV at /3 = 6.1, where the value of the 
string tension is obtained from the IP-IS mass 
splitting. Extrapolating linearly in to the con- 
tinuum limit, they obtain 



= 73 ± lOMeV. 



(4) 



This is still far short of the experimental value. 
However, it has been argued that the remain- 
ing discrepancy is consistent with the corrections 
arising from the quenched approximation S. 



0.00 



-0.02 




0.0 0.5 1.0 1.5 2.0 

1/Mp 

Figure 1. The hyperfine mass splitting My — Mp 
against 1/Mp, in lattice units, is shown both 
for heavy-heavy mesons (squares) , and for heavy- 
light mesons after extrapolation of ki to Kcrit = 
0.1431 (crosses). The fits to the heavy-heavy and 
heavy-light data are indicated by the solid and 
dashed lines respectively. 



Figure |i| shows My — Mp (where V and P rep- 
resent the vector and pseudoscalar respectively) 
for both heavy-light and heavy-heavy mesons ob- 
tained using the clover action on the 33 configu- 
rations. Fits are over the range 10 to 16, using 
the LS charm propagators. For both the heavy- 
light and heavy-heavy systems, linear behaviour 
with l/Mp is apparent; the non-zero intercept of 
the fits is presumably a manifestation of lattice 
artifacts. Setting = 2.73 GeV, we obtain 
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Table 1 

Masses (in lattice units) of the pseudoscalar and vector heavy-heavy mesons. 
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to be compared to the experimental values of 
146(1) MeV and 46(3) MeV respectively. An 
analysis using the clover action at [3 = 6.0, for 
Ki in the region of the strange quark mass, yields 
Tod* — niD = 76 1 21 @, clearly indicating a re- 
maining sensitivity to the lattice spacing, a. 



4. DECAY CONSTANT fp 

We extract the heavy-light pseudoscalar decay 
constant fp through a fit to the ratio. 



— tanhMp(T/2-<) (6) 
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where is the axial vector current, P — V'75^i 
and Z is defined through 

Z^{0\P\P). (7) 

S and L refer to quantities constructed from 
smeared and local quark propagators respectively. 
Zs and Mp are determined from the fit to 
{Psix,t)Pl{0)). The fits both to eq. |and to the 
pseudoscalar correlator are over the range 10 - 16. 
The quality of the data for eq. ^, together with 
the fits to the data, are illustrated in Figure |[ 
The expected behaviour from the HQET is 



fp^Mp^ A + 0{1/Mp) 
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where A is a constant, up to logarithmic terms. 
Figure H shows fp^/Mp against Mp for all ki and 
Kh, together with the value after the extrapola- 
tion of Ki to Kcrit = 0.1431 ai. 

We determine the behaviour of fpy/Mp with 
1/Mp, at Kcrit, from a correlated fit to 



fpy/Mp = a I 
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Figure 2. The ratio in eq. 6 is shown for all four 
Kh, at fixed ki = 0.14226. The solid lines repre- 
sent the fits to the data. 



In lattice units, a linear fit to the three heaviest 
points yields, 
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with x^/dof = 2.4/1, while a quadratic fit to all 
four points yields 

0.1201 § 



ai 

02 = -0.053 
ag = 0.0181 
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with x^/dof =0.1/1. A linear fit to all four points 
produces a noticeably poor fit, with x^/dof = 
12/2. 

Using the linear fit to the heaviest three points, 
we obtain 
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Figure 3. fpy^Mp against 1/Mp in lattice units 
is shown for all ki and Kh (squares). Also shown 
are the values after linear extrapolation of ki to 
Kcrit = 0.1431 (crosses). The solid and dotted 
lines represent a quadratic fit to all four points, 
and a linear fit to the heaviest three points re- 
spectively, at Kcrit- 
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GeV, 



while from the quadratic fit we obtain 
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where we take Za = 1 — 0.02^^, with g the bare 
coupling. The renormalisation coefhcient, Za, for 
the clover action is much closer to unity than that 
for the Wilson action {Zj ~ 1 - 0.132g2). Thus 
the uncertainty due to the choice of expansion 
parameter |11|, is only about 1% for the clover 
action, compared to about 10% for the standard 
Wilson action. 

The determination of fs is clearly more sensi- 
tive to the form of the extrapolation than that of 
/d. An analysis at /3 = 6.0 [|o|, using a hopping 
parameter expansion of the heavy quark propa- 
gator, favours a quadratic fit, yielding similar co- 
efficients (in physical units) to those quoted in 
eq. |ll|. The non-scaling corrections to Jb are ap- 



proximately 15% and 20% of the scaling expecta- 
tion, obtained from the constant term oi, for the 
linear and quadratic fits respectively. 

5. CONCLUSIONS 

The improvement in the measured value of 
the hyperfine splitting both for heavy-heavy and 
heavy- light mesons using an 0(a)-improved ac- 
tion is very encouraging, and the remaining dis- 
crepancies with experiment may be due to the 
effects of quenching. 

The preliminary results for and fs pre- 
sented here need further refinement, in partic- 
ular to take account of the uncertainties in the 
lattice spacing. Nevertheless, they indicate that 
non-scaling contributions to the HQET predic- 
tions are substantial, even at the mass of the B 
meson. It is particularly important to measure 
the static results, with light quarks computed us- 
ing the clover action, on the same set of con- 
figurations. A consistent extrapolation between 
the propagating quark results presented here and 
the static predictions would provide powerful evi- 
dence for the reliability of lattice studies of heavy 
quarks. 
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